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ABSTRACT 

We analyze more than 5000 days of high-quality Sun-as-a-star, radial velocity GOLF and photometric VIRGO/SPM helioseismic 
observations to extract precise estimates of the visibilities of the low-degree p modes and the ra-height ratios of the 1 = 2 and 3 
multiplets in the solar acoustic spectrum. The mode visibilities are shown to be larger during the GOLF red-wing configuration than 
during the blue-wing configuration, and to decrease as the wavelength of the VIRGO/SPM channels increases. We also show that the 
mode visibilities are independent of the solar activity cycle and remain constant overall with time, but that nevertheless they follow 
short-term fluctuations on a time scale of a few months. The 1=1 mode visibility also increases significantly toward the end of 
the year 1999. Comparisons with theoretical predictions are provided. Even though there is qualitative agreement, some significant 
discrepancies appear, especially for the 1 = 3 modes. The limb darkening alone cannot explain the relative visibilities of modes. These 
precise estimates should be used as references for the extraction of the p-mode parameters for any future investigation using the GOLF 
and VIRGO/SPM observations. 



the / = 2 and 3 multiplets (the so-called m-height ratio) dur- 
ing the peak-fitti ng procedure when e stimating the p-mode char- 
acteristics (e.g., ISalabert et al.L l2004l) . while the acoustic pow- 
ers of the / = 1,2, and 3 modes relative to the / = modes 
are allowed to vary (the so-called mode visibility). However, 
in asteroseismology, the mode visibilities are fixed to theoret- 
ical va lues dChristensen-Dalsgaard & Goughi 119821 : iPalle et all 
Il989ah because of the lower signal-to-noise ratio (SNR) of typi- 
cal data sets and shorter time series, and the m-height ratios are 
expressed as a function of the inclination of the rotation axis 
only (Appourchauxet al., 2008; Garcia et al., 200 9). Al though, 
the acoustic frequencies (e.g., Woodard & Novesj 1 1985b . pow- 
ers, and damping rates (e.g., IPalle et al.^ ,19901) vary with the 
magnetic activity cycle of the Sun and the first observation of 
similar variations of the p-mode parameters in a solar-like star 
was recently reported by .Garcia e t al. (2010), the mode visibil- 
ities are supposed to be independent of both solar and stellar 
magnetic activity. Frohlich et al.l (1 19971) observed variations with 
the height in the solar atmosphere in the intensity data collected 
by the Variability of Solar I Rradiance and Gravity Oscillation 
(VIRGO; Frohlich et all [1995 ) at the beginning of the Solar and 
Heliospheric Observatory (SoHO) mission. However, it remains 
to be verified if these values are different in the radial veloc- 
ity G lobal Oscillations at Low Frequency (GOLF; Gabrie let al.L 
'1995) measurements between the blue- and red-wing observing 
periods. 

The space-based Co nvection, R ot ation, and plane- 
tary Transits (CoRoT, 'Michel et^ l2008h and Kepler 
(Borucki et al., 2009) missions are currently providing in 
large quantities high-quality photometric measurements. After 
several years of Kepler observations of the same asteroseismic 
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1. Introduction 

The potential of helio- and asteroseismology to provide in- 
sights into the interiors of the Sun and other distant stars is 
unique. Measurements of the properties of the normal modes 
of solar oscillations have contributed grea tly to our knowl- 
edge of both the internal structure (e.g. Basu et all Il996t 
Antia & Chitre, 1998; Couvidat et al., 2003) and d ynamics (e.g . 
[Thompson et aUll996l;IChaplin et aUir999l:lGarcfa et al.Ll2004 
of the Su n's convective and radiative layers as far down to its 
core ( e.g. iMathur et all l2007t iGarcfa et all I2007L iMathur et all 
|2008|) . Moreover, strong constraints on the stellar properties 
have also been obtained with seismic analysis from both ground- 
based (e.g. iB azot et al, 2005) an d space-based (e.g. Brunti . 
120091; IChristensen-Dalsgaai-d et aTlllolotlMetcalfe et aUl2010l) 
\ observations. 

However, our knowledge of the properties of the solar and 
, stellar interiors depends on our ability to correctly measure the 
resonant modes of oscillations, which are commonly modeled 
as Lorentzian profiles in the acoustic power spectrum. The mode 
parameters (amplitude, linewidth, frequency, rotational splitting, 
and asymmetry) are then extracted by fitti ng the model to the 
data using maximum-likelihood functions (lAppourchaux et all 
Il998l) . The rotation of a star lifts the degeneracy of the os- 
cillation modes and splits the eigenfrequencies into 21 + I m- 
components, as -/ < m < +1, where I is the angular degree 
and m the azimuthal order. In Sun-as-a-star helioseismology, 
only the / + \m\ even components are visible and it is common 
practice to fix the height ratios between the m-components of 
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targets, the SNR will be high enough to measure the visibilities 
and m-height ratios in a wide range of solar-like stars in the 
Hertzsprung-R u ssell diagram at d i fferen t evolution stages 
dBedding et all 1201 Ot IChaplin et all l2010l) . The forthcoming 
observations from the radial velocity Stella r Observations 
Network Group (SONG, iGrundahl et"aU l2007l) and the photo- 
metric PLAnetary Transits and Oscillations of stars (PLATCQ) 
instruments will also provide extremely precise measurements 
of stellar oscillations in the very near future. Simultaneous 
radial velocity and photometric observations will be extremely 
useful in improving our understanding and modeling of stellar 
atmospheres. 

2. Data and analysis 

We analyzed simultaneous observations collected by the space- 
based, Sun-as-a-star Global Oscillations at Low Frequency 
(GOLF: lGabriel et all [1991 and Variability of Solar IRradiance 
and Gravity Oscillation ("VIRGO: iFrohlich et a l.. 1995) instru- 
ments onboard the Solar and Heliospheric Observatory (SoHO) 
spacecraft. The GOLF instrument is a resonant scattering spec- 
trophotometer measuring the Doppler wavelength shift - inte- 
grated over the solar surface - in the Di and D2 Fraunhofer 
sodium lines at 589.6 and 589.0 nm, respectively. The VIRGO 
instrument is composed of three Sun photometers (SPM) at 
402 nm (blue channel), 500 nm (green channel), and 862 nm 
(red channel). A total of 5021 days of radial velocity GOLF and 
intensity VIRGO/SPM observations starting on 1996 April 11 
and ending on 2010 January 8 were analyzed, with respective 
duty cycles of 95.4% and 9 4.7%. The GOLF v elocity time series 
were obtained following Garcia et al. ('2005) and calibrated as 
described in Jimenez-Reyes et al. (2003). Moreover, in the con- 
text of this paper, it is important to remember that the GOLF 
signa l has a small intensity pollution of around 14% (iPalle et al.L 
. For technical reasons, GOLF has been observing only one 
side of the sodium doublet: in the blue wing from April 1 1 1996 
until June 25 1998 and again from November 18 2002 until now: 
and, in the meanwhile, between September 1998 until November 
1 8 2002, the measurements were obtained from the red wing of 
the sodiu m doublet, wh i ch o riginates higher up in the solar at- 
mosphere (iGarcfa et al.L 1200 5). In addition, two extended gaps 
are present in the analyzed time series because of the temporary 
loss of the SoHO spacecraft. The first gap of ~100 days hap- 
pened during the summer of 1998 and was due to a failure in the 
gyroscopes. The second one of a period of ~1 month in January 
1999 occurred while new software was being uploaded to the 
spacecraft. 

The power spectrum of each time series was fitted to esti- 
mate the mode parameters of the / - 0, 1,2, and 3 modes using 
a multi-step iterative method (ISalabert et al.L l2007l) . Each mode 
component of radial order n and angular degree / with azimuthal 
ord er m was parameterized with an asymmetric Lorentzian pro- 
file (iNigam & Kosovichevi[T998h . as 

(1 +ai^nXl,m,nf + Q"^,, 
-Cl,m,iM - Hi,„„ ; '—, (1) 

1 + Xf 

t,m,n 

where xi^„^„ = 2(v - v/,„,„)/r,,„, and v,,,„,„, F,,,,, and rep- 
resent the mode frequency, the linewidth, and the height of the 
spectral density, respectively. The peak asymmetry is described 
by the parameter q;;,„. The rotation is treated using a first-order 
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approximation (iLedouxi Il95lh . thus the mode frequency is ex- 
pressed as v; ,„_„ = v/,„ - w5v/ „, where v/ „ is the central frequency 
and 6vi „ is the rotational splitting. We note that we assumed the 
same linewidth for the components of a given multiplet. Because 
of their close proximity in frequency, modes are fitted in pairs 
(i.e., I = 2 with 0, and / = 3 with 1). While each mode pa- 
rameter within a pair of modes are free, the peak asymmetry is 
set to be the same within pairs of modes. When visible in the 
power spectra, the I = 4 and 5 modes were included in the fitted 
model. An additive value B is added to the fitted profile repre- 
senting a constant background noise in the fitted window. Since 
SoHO observes the Sun equatorwards, only the l + \m\ even com- 
ponents are visible in Sun-as-a-star observations of GOLF and 
VIRGO/SPM. To reduce the number of free parameters and sta- 
bilize the peak-fitting procedure, it is common practice to fix 
the height ratios between the visible m-components for both the 
/ = 2 and / = 3 multiplets to their estimated values. To derive 
observational estimates of the m-height ratios, the m = +2 and 
m = components of the I = 2 multiplets, and the m - +3 
and m = +1 components of the / - 3 multiplets were fit- 
ted using independent heights Hi,„^„, assuming that components 
with opposite azimuthal order m have the same heights, i.e., 
Hi-m,n - Hi^+ii,^„. Finally, the mode parameters were extracted 
by maximizing the likelihood function by assuming the power 
spectrum statistics is given by ax^ with two degrees of freedom. 
The natural logarithms of the mode height, linewidth, and back- 
ground noise were varied resulting in normal distributions. The 
formal uncertainties in each parameter were then derived from 
the inverse Hessian matrix. The blue and red periods of GOLF 
were also analyzed separately, as well as the mean power spec- 
trum of the three VIRGO/SPMs (blue, green, and red channels). 

3. Mode visibilities 

Because of averaging effects over the solar disk with Sun-as-a- 
star observations, the observed relative amplitudes of the modes 
are modified by visibility factors V/. It is possible to recover the 
visibilities by assuming that the intrinsic amplitudes of modes 
with close frequencies are similar. By denoting A/„ to be the 
observed amplitude of a mode, we find that Vi/Vo ~ A/_„/A/=o,«' 
when the modes (/, «) and (/ = 0, n') have close frequencies. 
Thus, we can measure the visibilities in the spectrum as 

(Vi/Vof ^Pi.JPo,„, (2) 

iV2/Vof ^P2,,-l/P0,n, (3) 

and 

iVi/Vof ^P3,„-l/Po,n, (4) 

where P/ „ is the mode power, directly linked to the square of the 
mode amplitude, defined as P/ „ oc //; „r/ „, while the height Hi „ 
of a given multiplet (I, n) is defined as the sum of the heights of 
its m-components, as: 

Hl,„ - ^ Hijn,n- (5) 

m=-l 

Since F/,,, is assumed to be a smoothly varying function also 
of the frequency, we could directly use the ratios of the mode 
heights Hi„ to measure the visibility. We decided to use the 
mode powers instead because, as demonstrated by Bal lot et al.l 
(120081) . they can be accurately measured even when the individ- 
ual m-components within a given multiplet start blending around 



2 



D. Salabert et al.: Mode visibilities in Sun-as-a-star helioseismic observations 



2 . 5 



1 . 5 



0.5 




1500 2000 2500 3000 3500 

Frequency (|IHz) 



3r 

2.5- 

2 

1 . 5 

1 

0.5 



4000 



Corrected 




1500 2000 2500 3000 3500 

Frequency (flHz) 



4000 



Fig. 1. Raw (left) and corrected (right) mode visibilities (Vi/Vo)^ of Z = 1 (blue dots), / - 2 (red squares), and I - 3 (black diamonds) 
relative to Z = as a function of frequency in the case of the radial velocity GOLF (solid lines) and intensity VIRGO/SPM (dashed 
lines) observations. 



3100 yuHz because of their increasing linewidths with frequency. 
Although consistent values are obtained between the two defini- 
tions, the computation of the visibilities defined as the ratios of 
the mode hei ghts is hmited t o frequ encies lower than 3100 fjHz, 
as shown bv ISalabert et al.l (1201 1 ^. The / = 1,2, and 3 mode 
visibilities calculated as in Eqs. |2] [3] and |4] are represented in 
Fig. [T]as a function of frequency, for both radial velocity GOLF 
and intensity VIRGO/SPM (i.e., the mean power spectrum of 
the three SPMs on VIRGO) observations. These raw mode vis- 
ibilities (left panel of Fig. [U vary with frequency - especially 
for the I - 1 mode - that is mainly due to the variation in the 
mode heights and linewidths with frequency. The heights and 
linewidths vary sensitively with frequency over half the large 
frequency separation (~ 67yuHz), i.e., the separation in frequency 
between the I = and I - 1 modes of the same radial order n. 
Thus, the visibilities are biased and to take account of this ef- 
fect, this frequency dependence was corrected by interpolating 
with a spline function the powers f/ „ of the I = 1, 2, and 3 
modes to the frequencies of the / = mode. The corrected mode 
visibilities are represented in the right panel of Fig. \T\ As the 
corrected mode visibilities calculated using the mode powers do 
not display significant variations with frequency, we can safely 
average them over frequency to obtain estimates of the true val- 
ues. We note that we use the expression mode visibility here to 
refer to the corrected mode visibility. Figure |2] shows the mode 
visibilities (Vi/Vof, (Vi/Vof, and (Vs/Vof averaged over fre- 
quency, using the errors as weights, measured from GOLF and 
VIRGO/SPM observations as a function of the angular degree 
/. The represented values are tabulated in Table [1] The GOLF 
visibilities were averaged between 1800 and 3600 yuHz for all 
/ values, and the VIRGO/SPM visibihties from 2200 to 3600 
fiHz for the I - 1 and 2 modes (and from 2500 yuHz for the 
1-3 mode). The lower limits are defined by the SNR of the 
modes, the intensity measurements such as VIRGO/SPM hav- 
ing a higher noise level at low frequency and less sensitivity to 
the modes. The mode visibilities in the blue-wing and the red- 
wing GOLF configurations, as well as the visibilities measured 
separately in the blue, green, and red VIRGO/SPM channels 
are also given in Table [1] Significant differences are observed 
depending on the observing wing configuration of GOLF and 
the VIRGO/SPM channels analyzed. The mode visibilities are 
greater during the GOLF red-wing configuration than during the 




1 2 
Angular degree 1 

Fig. 2. Mode visibilities (Vi/Vq)^ averaged over frequency as a 
function of angular degree I for the radial velocity GOLF (black 
diamonds, solid line) and intensity VIRGO/SPM (red squares, 
dashed line) measurements. 

blue-wing configuration, the / = 3 mode being for instance twice 
as visible during the red-wing period than in the blue-wing pe- 
riod. The mode visibilities are shown to decrease as the wave- 
length of the VIRGO/SPM channels increases. 

4. Height ratios between m-components 

The height ratios /?/ „ between the visible m-components of the 
I = 2 and / - 3 multiplets in Sun-as-a-star observations are 
defined, respectively, as 



^2,n - Hi=2,m=0,n / Hi=2,m=±2,n, 

and 

^3,« - Hi=2,m=±l,n/Hl=3,m=±3,n- 



(6) 



(7) 



As the linewidths are forced to be the same between m- 
components of a given multiplet (I, n) in the peak-fitting proce- 
dure, the ratios thus measured using the heights are iden- 
tical to the ones calculated using the powers P/_,„,„. These ratios 
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Table 1. Mode visibilities in the radial velocity GOLF and intensity VIRGO/SPM measurements. 



(V,/Vo)' 


GOLF 


GOLF 
Blue wing 


GOLF 
Red wing 


VIRGO/SPM 


VIRGO/SPM 
Blue 


VIRGO/SPM 
Green 


VIRGO/SPM 
Red 


(V2/V0)' 


1.730 ±0.017 
0.884 ±0.012 
0.167 ±0.002 


1.707 ± 0.023 
0.827 ± 0.017 
0.138 ±0.003 


1.752 ± 0.028 
1.016 ± 0.022 
0.259 ± 0.006 


1.547 ±0.016 
0.655 ±0.012 
0.109 ±0.003 


1.578 ± 0.016 
0.696 ± 0.012 
0.122 ± 0.003 


1.540 ±0.016 
0.632 ± 0.011 
0.094 ± 0.002 


1.354 ±0.014 
0.465 ± 0.008 
0.063 ± 0.002 



Notes. The mode visibilities in the blue-wing and the red-wing GOLF configurations are also given, as well as the visibilities measured separately 
in the blue, green, and red VIRGO/SPM channels. 
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Table 3. Modeled visibilities and m-height ratios computed for 
the three VIRGO/SPM channels. 
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Fig. 3. m-height ratios Pi of the / = 2 (top) and I - 3 modes 
(bottom) as a function of frequency extracted from the GOLF 
radial velocity measurements. 

Table 2. Height ratios fii between the m-components of the I - 
2 and / = 3 modes measured in the radial velocity GOLF and 
intensity VIRGO/SPM observations. 



Pi 


GOLF 


VIRGO/SPM 


Pi 


0.634 ± 0.033 


0.751 ±0.059 


P3 


0.400 ± 0.020 


0.633 ± 0.065 



Pi^„ are represented in Fig. |3] as a function of frequency in the 
case of the GOLF observations. Some marginal variations with 
frequency can be seen in Fig. |3] but we assume that the height 
ratios are independent of frequency, and that an average over the 
spectrum will yield an estimate of the true value. The m-height 
ratios of the I = 2 and / = 3 multiplets averaged over frequency, 
using the errors as weights, are given in Table |2] for GOLF and 
VIRGO/SPM observations. No significant differences are found 
in the m-height ratios between the GOLF blue and red observing 
wings, and also between the three VIRGO/SPM channels. 

5. Modeling visibilities 

5.1. Photometric observations 

As previously mentioned, the apparent amplitude value of a 
mode depends on the way that the contributions of each element 
of the solar disk to the total signal are averaged. When the contri- 
bution of a solar-disk element to the total flux depends only on its 
distance to the limb, the mode visibility Vi and the m-height rati o 
are decoupled re.g.. lGizon & SolankiL[2003HBallot et aUl2006h . 
The observed amplitude A^^°^2 of a mode (n, /, m) is linked to the 
intrinsic amplitude A„ by means of the relation 



Mode visibility & 


Blue 


Green 


Red 


m-height ratio 


(403;um) 


(501yum) 


(863;um) 


(Vi/VoY 


1.60 


1.54 


1.45 




0.66 


0.57 


0.45 


(V3/V0)' 


0.06 


0.03 


0.01 


P2 




0.667 




Pi 




0.600 





Notes. The m-height ratios do not depend on the channel. 



where 



V, 



and 



y/2n{2l+l) f WijiWM^d^, 
Jo 



{l + \m\)V ' 



(9) 



(10) 



where Vi denote the Legendre polynomials, !P"' are the associ- 
ated Legendre functions, and / is the inclination of the rotation 
axis relative to the line of sight, which is / - 90° for the Sun 
to a very good approximation. The function W(/i) is the weight- 
ing function, which depends only on the distance to the limb yu 
classically defined as the cosine of the angle MCO, where M is 
the considered point at the solar surface, C the center of the Sun, 
and O the observer. Within this approximation, and by assum- 
ing that there is an equipartition of the energy between modes 
of close frequency and that the mode width depends only on the 
frequency, we find that 



//,///„ = PiIPq = {ViIVqY, and 

Pl,m'/Pl,m - Hum' li^lm - {ai,m' 1 0!lmf' . 



(11) 

(12) 



n.Lin 



ViA„,,, 



(8) 



For photometric measurements such as VIRGO/SPM observa- 
tions, the initial assumption - W depends on fi only - is veri- 
fied because the contribution depends mainly on the limb dark- 
ening. As a consequence, the amplitude ratios between the m- 
components are purely geometrical effects, and the mode visi- 
bility Vi is computed by assuming the weighting function to be 
W(/i) - h{iJ)IUQ), the limb-darkening profile at the observed 
wavelength A. We computed the expected visibility by consid- 
ering the spectral bands of the three VIRGO/SPM channels . 
We used the limb-darkening law from iNeckel &Labl ([T99l . 
Results are given in Table |3] 

We compare these predictions to observations (cf. Table [TJ. 
First, the variations with color are correctly reproduced for the 
different modes. The computations recover the observations to 
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within a few percent accuracy, except for the / - 3 modes, which 
are systematically higher than expected. Nevertheless, we note 
that, even for I = 1 and 2 modes, half of the predictions have 
discrepancies with the observations that are larger than 3 cr. In 
addition, the error bars in Table [1] are purely statistical, and do 
not take into account methodological or systematic errors. The 
error bars might then be underestimated. There are several pos- 
sible explanations of these differences. First, the systematically 
larger amplitudes observed for the / = 3 modes could be due 
to the intrinsic amplitudes being larger than expected. It can- 
not be excluded, but is hardly justifiable from a theoretical point 
of view today. Secondly, the differences might come from the 
modeled weighting function being unable to perfectly reproduce 
the real one. The limb-darkening law coming directly from ob- 
servations is reliable. By identifying the limb-darkening profile 
with the weighting function, we assume that the relative inten- 
sity fluctuations 61 /I induced by a mode are the same every- 
where on the visible solar disk. However, the limb is darker be- 
cause the light comes mainly from higher layers of the photo- 
sphere. With an adiabatic approximation, and knowing that the 
temperature T varies slowly in the atmosphere, it is shown that 
= 46T/T oc where is the displacement. Since the 
mode is evanescent in the atmosphere. 



Table 4. Modeled visibilities and m-height ratios in radial ve- 
locity for GOLF observations. 



P^r exp 



•(l -LLt/colf'^ h/H 



(13) 



where h is the height in the atmosphere, H is the pressure 
scale height, and a>r is the cu t-ofi" frequency (see for instance 
IChristensen-Dalsgaard L l2003l) . As the photosphere is a thin 
layer, we generally assume that does not vary too much. 
Nevertheless, because of the strong drop in density p, it is pos- 
sibly not exactly true. Moreover, non-adiabitic effects are proba- 
bly not negligible. Since Vi=3 modes appear to be relatively more 
sensitive to the shape of Wi/J.) close to the limb, the observed dis- 
crepencies could be due to the approximated physics of modes 
in the atmosphere. 

Concerning the m-height ratio, the value are close to the ob- 
served ones (cf. Table |2]i, even if, for 1 = 2 modes, the observed 
central component is marginally larger than predicted (~ 1 Act). 
Since this term is normally insensitive to the physics of the solar 
atmosphere, discrepencies can be due to differences in the intrin- 
sic amplitudes of the axisymmetric component. Nevertheless, 
GOLF does not have similar systematics (see next section), 
meaning that these diff'erences are not conclusive. 

5.2. GOLF observations 

For GOLF observations, the weighting function does not depend 
only on yu, and we have to perform complete computations tak- 
ing into account the instrumental response, which depends on 
the positions in the NaDi and D2 lines that are observed. Thus, 
mode visibilities change according to the observed wing (blue 
or red), but also with the epoch of the observations: according to 
the position in the SoHO orbit, the Doppler shift change and 
the lines are probed at different heights in the w ings. Details 
of the instrumental response can be found in either' Garcia et alJ 
( Il998l) or Ulrich et al. (2000). We numerically computed the in- 
tegrals over the whole disk for the different modes through a 
simulated instrumental response taking into account a simplified 
limb-darkening of the NaD lines, the solar rotation and differ- 
ential rotation, and the gravitational reddening. We first assume 
that the modes induce purely vertical motions. 

The results of these computations are listed in TableU] First, 
we note that the m-height ratio does not depend noticeably on 



Mode visibility & 


GOLF 


GOLF 


m-height ratio 


Blue wing 


Red wing 


(Vi/VoY 


L84 


L86 


(V2/V0)- 


LOS 


L13 




0.27 


0.30 


Pi 


0.59 


0.58 


Pi 


0.43 


0.40 



the considered wings and is in good agreement with the observed 
ratio (see Table |2]i. For the visibility, the global trend is recov- 
ered, but the difference between blue and red wings for GOLF 
is larger than expected: the observed value for the blue wing is 
significantly smaller than expected. 

We also performed complementary computations to take 
into account the horizontal motions of the modes. Since 
the contribution of the horizontal motions increases when 
the frequency decreases and the mode degree increases (e.g. 
Christensen-Dalsgaard , l2003h . it generates a small variation 
in the mode visibility with frequency by a few percent. 
Nevertheless, it is negligible compared to the observation errors. 



6. Temporal evolution of the mode visibilities 

As for the p-mode frequencies (e . g., Woodard & Novesl 1985L 
iPalle et all [19895 IChaplin et al.L 120071: ISalabert et al.L l2009^ ■ 
the acoustic power and damping parameters of the low-degree 
acoustic rn o des ha v e been found t o vary with solar activi t y (e.g. , 
Palle et all IT990t IChapUn et all I2OO0I ISalabert et all l2003l 
Simoniel lo et al However, it has never been demon- 

strated whether the mode visibilities depend or not on the Sun's 
magnetic activity. To verify this, we analyzed the GOLF and 
VIRGO/SPM observations by dividing the original time series 
into contiguous 365-day and 9L25-day non-independent sub- 
series (shifted by 1/4) and fitting the power spectrum of each 
subseries. For the SoHO vacation in 1998, the subseries with 
duty cycles less than 50% were removed from the analysis. 
The mode parameters were extracted in the same manner as ex- 
plained in Sec. |2] with the only difference that the height ratios 
between the visible m-components of the I = 2 and 3 modes were 
fixed to the values determined in Sec. |4] and given in Table |2l in 
order to stabilize the fitting procedure of these short time series. 
The mode visibilities obtained as in Eqs. |2l [3] and |4] for each 
power spectrum were corrected for their frequency dependences 
by interpolating the powers of the / = 1,2, and 3 modes to the 
I = frequencies as described in Sec. [3] and their weighted av- 
erages over frequency obtained. 



6. 1 . The 365-day subseries 

Figure|4]shows the temporal evolution of the / = 1,2, and 3 mode 
visibilities measured in the non-independent 365-day GOLF and 
VIRGO/SPM power spectra. The mode visibilities do not vary 
with solar activity, and remain constant overall over the years, 
within the uncertainties. However, a clear dependence on the 
GOLF wing configurations, as indicated by the vertical dashed 
lines, can be discerned, mainly in the / = 3 mode. The GOLF vis- 
ibilities increase during the red-wing observing period between 
the years ~ 1998 and ~ 2002, which is consistent with the re- 
sults in Table[T] On the other hand, the VIRGO/SPM visibilities 
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do not exhibit these changes. However, the VIRGO/SPM visibil- 
ity of the / - 1 mode has a relatively large increase centered on 
the year 2000, which is also present in the GOLF data in addi- 
tion to the larger visibilities during the red-wing operation. The 
VIRGO/SPM I - 2 mode visibility increases over the same pe- 
riod of time appearing to follow the GOLF visibility. As the level 
of precision of the data, the VIRGO/SPM / = 3 visibilities do not 
show any changes comparable to the ones observed in / = 1 and 
2. The weighted mean values over time of the GOLF visibilities 
are 1.747 + 0.011, 0.862 + 0.008, and 0.171 + 0.003 for / = 1, 
2, and 3 respectively. In the case of VIRGO/SPM, the values are 
1.605 + 0.012, 0.641 + 0.008, and 0.089 + 0.004, respectively. 
These results are in agreement with Table[T] Similar conclusions 
are obtained with the individual VIRGO/SPM channels. 

6.2. The 91 .25-day subseries 

Figure |5] shows the temporal evolution of the / - 1 and 2 mode 
visibilities measured in the non-independent 91.25-day GOLF 
and VIRGO/SPM power spectra. We do not show the results for 
the I = 3 mode as they have a much larger scatter in as short 
time series as these, because of the mode's smaller visibility. 
Although, Fig.|5]confirms that the mode visibilities are indepen- 
dent of the solar activity, it shows fluctuations over short peri- 
ods of time, of about 6 months. The same analysis performed 
using different lengths of subseries and amounts of shift (1/2, 
1/3, and 1/4) confirms that this six-month period is not an arti- 
fact of the methodology. Moreover, the analysis of the artificial 
SolarFLAG data (Chaplin et al., 2006) does not detect this peri- 
odicity. However, its origin remains unclear. Some of these fluc- 
tuations have significant increases, common to the GOLF and 
VIRGO/SPM measurements. The enhancement of the I - 1 vis- 
ibility observed in the 365-day subseries could be due to a se- 
ries of higher-than-average peaks between 1999 and 2001 that 
are observed in the 91.25-day visibilities. Inspection of the in- 
dividual mode powers seems to indicate a significant reduc- 
tion in the / = power during that period. Similar results are 
obtained with the individual VIRGO/SPM channels. Thus, to 
check whether the origin of this increase in the mode visibil- 
ities is genuine and unrelated to intrinsic changes in SoHO, 
we performed in Appendix |A] a similar analysis using the in- 
dependent observations collected by the ground-based, multi- 
site network B i rming ham Solar Oscillations Network (BiSON; 
lElsworth et al.L 1995h and G lobal Oscillation Network Group 
(GONG: lHarvev et al Lll996l) . We show in Appendix|B]that iden- 
tical temporal fluctuations seen in Figs. |4]and|5] including the 
increase in the / - 1 visibility in 2000, are also present in the 
BiSON and GONG observations. 



7. Conclusions 

We have analyzed more than 5000 days of high-quality helio- 
seismic observations collected by the space-based, Sun-as-a-star, 
radial velocity GOLF and intensity VIRGO/SPM instruments 
onboard the SoHO spacecraft to extract precise measurements of 
the low-degree, p-mode visibilities and m-height ratios. The vis- 
ibilities of the / = 1, 2, and 3 modes relative to the I - mode are 
defined as the ratios P/=i/P;=(), Pi=2/Pi=o, and P/=3/P;=o respec- 
tively, where Pi is the acoustic power of a given mode of angu- 
lar degree I. The height ratios between the m-components of the 
1 = 2 and / = 3 multiplets are defined as ///=2,m=o/^/=2,m=±2 and 
Hi=i,m=±\ I Hi=3,m=±i respectively, where Hi ,„ is the height of the 
individual m-components of the mode. We derived observational 



values that we then compared to theoretical values computed us- 
ing the instrumental responses of GOLF and VIRGO/SPM by 
taking into account limb-darkening functions. In general, the 
predicted values qualitatively agree with the observations. We 
also showed that the mode visibility is a function of the wing 
configuration of GOLF, the / = 1,2, and 3 modes being more vis- 
ible during the red-wing period than during the blue- wing one. 
The mode visibilities were also shown to depend on the wave- 
length of the VIRGO/SPMs, as predicted. However, some dis- 
crepancies between observations and theory remain, especially 
for the I = 3 modes. These differences imply that the limb dark- 
ening alone is insufficient to explain the observed visibilities. 
This has to be analyzed in depth to be able to make reliable 
predictions for other stars. Nevertheless, we note that the dis- 
crepencies identified here are smaller than the error bars in mode 
amplitudes obtained with a few months of CoRoT or Kepler ob- 
servations. 

By analyzing 365-day and 91.25-day GOLF and 
VIRGO/SPM subseries, we also demonstrated for the first 
time that the mode visibilities are constant in time and do 
not vary with the solar magnetic activity. However, the mode 
visibilities exhibit short-term temporal variations of about 
six months that are common to the GOLF and VIRGO/SPM 
observations. Some of these fluctuations increase significantly, 
and a significant increase in the / = 1 mode visibility was found 
at the end of the year 1999. The origin of both these short-period 
fluctuations and this bump at the end of 1999 is unclear but 
might be related to the stochastic excitation of the individual 
modes or the variations in the acoustic power and damping rate 
of the solar oscillations with time. 

In the future, the observations from the radial velocity SONG 
and photometric PLATO instruments will provide invaluable in- 
puts to model the stellar atmospheres. Finally, these newly con- 
strained / = 1,2, and 3 mode visibilities and m-height ratios 
provided in this paper should be used as inputs to any investi- 
gation of the Sun using the radial velocity GOLF and intensity 
VIRGO/SPM observations. 
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Appendix A: Mode visibilities in BiSON and GONG 
integrated observations 

This analysis presented in these appendices was motivated by 
the observation of an increase in the temporal evolution of the 
visibilities, mainly for the I - 1 mode, at the end of the year 
1999, seen in the Sun-as-a-star GOLF and VIRGO/SPM in- 
struments onboard the SoHO spacecraft. We wish here and in 
Appendix iBl to establish whether the origin of these variations 
might come from SoHO itself and not from the Sun. We de- 
cided then to perform the same analysis using two indepen- 
dent helioseismic datasets collected by the ground-based, multi- 



Table A.l. Mode visibilities extracted from the analysis of the 
BiSON and integrated GONG time series. 



(Vi/VoY 


BiSON 


GONG 


(Vi/VoY 
(Vi/VoY 


1.895 ±0.016 
1.296 ± 0.016 
0.467 ± 0.006 


1.606 + 0.014 
0.652 ± 0.008 
0.028 ± 0.001 



site BiSON, and GONG networks. The unresolved BiSON and 
spatially-resolved GONG networks are both composed of six 
stations located at selected longitudes around the world. The 
instruments at each BiSOl^ site make Sun-as-a-star observa- 
tions of the Doppler shi ft of t he potassium Fraunhofer line at 
770 nm ( Elsworth et all [T995I) . while the GONG cameras use 
Michelson Doppler interferometer-base d instruments measu ring 
in the absorption line Ni I at 676.8 nm (iHarvev et all Il996ll . In 
these appendices, we used the spatially-integrated GONG time 
series, which is analogous to the unresolved observations of 
GOLF, VIRGO/SPM, and BiSON, although the visibilities for 
the higher degree modes are different, with the relative strength 
falling off faster as a function of I in the GONG data. A total 
of 5339 days, starting on 1995 January 1, of BiSON observa- 
tions, and of 5112 days, starting on 1995 May 7, of the GONG 
integrated data, with respective duty cycles of 8 1 .0% and 85.4%, 
were analyzed in the same manner as described in Sec. [2] The 
daily temporal sidebands due to the diurnal gaps in ground-based 
observations were included in the fitting model (Sec. |2]l. The 
BiSON and GONG integrated I - 1,2, and 3 mode visibiHties, 
obtained as explained in Sec.[3j are represented as a function of 
frequency on Fig. lA.ll and given in Table lA.ll once averaged 
over frequency between 1800 and 3100 yuHz. The / - 3 mode 
has a very small visibility in the GONG integrated data and was 
averaged from 2500 yuHz only, due to its lower SNR at low fre- 
quency. The ;7i-height ratios were also estimated and are given 
in Table lA.2] Chaplin et al. ( 2001) calculated the w-height ratios 
in the BiSON data finding 0.55 + 0.04 and 0.38 + 0.02 for the 
/ = 2 and / = 3 modes, respectively. Our analysis returns larger 
ratios than theirs, but are consistent within Icr for I = 2 and 2cr 
for 1 = 3. These differences might be cause by differences in the 
total length of the datasets used, 8 years in Chaplin et al. (200 ij) 
compared to 14 years here, meaning a higher frequency resolu- 
tion in the present case and different periods of time analyzed - 
indeed, variations with solar activity cannot be totally excluded. 
We calculated these height ratios over the same period of time 
as in ChapHn et al. (2001) and found 0.57 + 0.04 and 0.4 ±0.03 
for / - 2 and 3, respectively, which agree with iChaplin et al.l 
(1200? ) to 1 cr. The remaining differences might be attributable to 
the peak-fitting itself, such as the size of the fitting windows, or 
the small multiplet frequency asymmetry at Z = 2 and 3, which 
is not taken into account in our analysis. In the case of the in- 
tegrated GONG time series, no previously measured m-height 
ratios were found in the published litterature. We note that since 
the visibility of the 1 = 3 mode is close to being zero in the 
GONG integrated data, the m-height ratio for I = 3 could not be 
measured accurately. 



- Respective datasets are available from 

http : //bi son ■ ph ■ bham . ac . uk/data . php | and 
http : / /gong . nso . edu/ data/. 
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Fig. A.l. Mode visibilities of I - 1 (blue dots), I = 2 (red squares), and I - 3 (black diamonds) relative to Z = as a function of 
frequency extracted from the analysis of the BiSON (left) and GONG integrated (right) observations. 
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Fig. A.l. Temporal variations in the I - 1 mode visibilities extracted from the analysis of the non-independent 365-day (left column) 
and 91.25-day (right column) power spectra of the GOLF, VIRGO/SPM, GONG, and BiSON (from top to bottom) instruments. The 
vertical dashed lines correspond to 1999 October 10. 



Appendix B: Temporal evolution of the / = 1 mode 
visibility in several Sun-as-a-star helioseismic 
instruments 

We compare the temporal evolution of the low-degree mode vis- 
ibilities of the solar oscillations measured with the spaced-based 



9 



D. Salabert et al.: Mode visibilities in Sun-as-a-star helioseismic observations 

Table A.2. Height ratios /?/ between the m-components of the 
/ - 2 and I = 3 modes measured in the BiSON and integrated 
GONG time series. 





BiSON 


GONG 




0.625 ± 0.032 


0.882 + 0.045 




0.445 ± 0.020 





Notes. Owing to the very small visibility of the / = 3 mode in 
the GONG integrated data, the m-height ratio for / = 3 could not be 
accurately measured. 

GOLF and VIRGO/SPM instruments onboard the SoHO space- 
craft, to the visibilities measured with the ground-based, multi- 
site BiSON, and GONG networks. The BiSON and GONG in- 
tegrated time series were divided into contiguous 365-day and 
91.25-day non-independent subseries (shifted by 1/4) and ana- 
lyzed in the same manner as explained in Sec. |6] Figure IA.2I 
shows the temporal evolution of the visibility of the I - 1 mode 
extracted from the analysis of the non-independent 365-day and 
91.25-day GOLF, VIRGO/SPM, BiSON, and GONG integrated 
power spectra. The bump of the / - 1 mode visibility at the end 
of the year 1999 observed in the 365-day subseries (see Fig. |4]i 
is present in all the analyzed independent datasets. As already 
mentioned in Sec. |6l this bump is most likely the result of a 
series of higher-than-average peaks between 1999 and 2001 ob- 
served in the 91.25-day visibilities. Moreover, similar short-term 
variations of about 6 months can be observed in the 91.25-days 
subseries of the GOLF VIRGO/SPM, BiSON, and GONG in- 
struments, some of them exhibiting significant and common in- 
creases. The same analysis performed using different lengths of 
subseries and overlapping factors confirms that this six-month 
periodicity is not an artifact of the methodology. 



10 



